Thermal quenching of electronic shells and channel competition in cluster fission 
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Experimental and theoretical studies ol fission ol doubly- 
charged Li, Na, and K clusters in the low fissility regime reveal 
the strong influence ol electronic shell elFects on the flssion 
products. The electronic entropy controls the quenching ol 
the shell effects and the competition between magic-lragment 
channels, leading to a transition Irom lavored channels ol 
higher mass symmetry to the asymmetric channel involving 
the trimer cation at elevated temperatures. 
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Explorations of factors controlling the stability and the 
collapse or fragmentation of finite systems range from 
early studies of hydrodynamical instabilities and droplet 
formation [1] to studies of fission and fragmentation of 
atomic nuclei [2] , and of atomic and molecular free clus- 
ters and supported nanostructures [3,4]. Here we focus 
on the charge instability and fission of doubly charged 
cationic metal clusters, that is — > Mp -I- M^_p. 
Such instabilities in macroscopic fluid systems have been 
studied by Rayleigh [1], and adaptations of Rayleigh's 
model (RM) for the description and analysis of nuclear [2] 
and cluster fission [5] are referred to as the Liquid Drop 
Model (LDM). hi the LDM, the clusters are viewed as 
charged classical liquid drops whose shape and dynamics 
are controlled by the competition between the repulsive 
Coulomb (-Bcouiomb) and the binding surface (i^surfaco) 
energies [6]. 

Various fission regimes may be classified by the fissil- 
ity parameter X = £^couiomb/2^^surfacc';^here the super- 
script "sph" indicates that the equilibrium shapes of the 
LDM are spherical. For X > 1, spontaneous barrier- 
less fission is predicted and observed [7,8]; this is also 
the regime studied originally by Rayleigh who predicted 
a mass-symmetric fission mode (see also Ref. [7]). In 
the X < \ regime, which is the focus of our paper, the 
fissioning system must overcome a barrier. For atomic 
clusters in this regime, tunneling is suppressed, unlike 
the case of nuclei, and fission requires thermal activation. 
From Xcr — 1, the RM critical size for barrierless fission 
is = z^e^/(167rrgCr), where Ts is the Wigner-Seitz 

radius and cr is the surface-energy; for doubly charged 
(z = 2) Li, Na, and K clusters, one has Nf^ « 10. 
Thus, with the exception of the smallest ones, fissioning 
clusters must be hot [9] in order to overcome the fission 



barrier. 

For the low fissility X < 1 regime, the LDM predicts 
the dominance of a strongly mass asymmetric fission 
process, as deduced from the finite-temperature [6,10] 
Qp = Fp + F^_p — F"^ values, that express the free- 
energy balance for the possible (P, N — P) fission chan- 
nels. Such asymmetric fission has been observed in exper- 
iments for Na^ with iV 24 (i.e., near the appearance 
size, see below), where the preferred fission channel in- 
volved the "magic" Naj^ fragment [9]. The prevalence 
of this channel, compared to other asymmetric ones in 
the neighborhood of the LDM minimum, originates from 
quantum mechanical shell effects [11]. 

In this paper we present experimental and theoretical 
results, pertaining to cluster fission for X << 1, which 
cannot be described by the classical RM, demonstrat- 
ing instead the dominance of quantum-size effects, i.e., 
electronic shell contributions [6]. In particular, we study 
alkali-metal (Li, Na, and K) clusters with N = 24, 26, 28, 
and 30. These systems correspond to the lowest observ- 
able fissility regime with N » N^^; e.g., for iV = 30 
one has X « 0.25 for the doubly charged alkali clusters. 
Moreover, this range of sizes is particularly suitable for 
explorations of shell effects, since it includes cases where 
the clusters may fission into doubly "magic" fragments, 
i.e., Ml+ M+ + M+ and M+ + M+. 

We find that, in addition to M;^, the more mass sym- 
metric magic fission channels (involving the magic frag- 
ments M(j" and M^j^) do appear and compete with the 
channel involving the trimer (even becoming the pre- 
ferred channels in several instances). Furthermore, from 
comparisons between the experiments and calculations 
based on the finite-temperature shell-correction method 
(FT-SCM) [12], we find that the channel distributions 
are greatly influenced by thermal quenching of the shell 
effects due to the increase (compared to T = 0) of the 
electronic entropy; this quenching results from thermal 
promotion of electrons to unoccupied single-particle lev- 
els which smears out the shell structure. The more mass- 
symmetric fission channels are favored at lower temper- 
atures, whereas the asymmetric ones (e.g., involving the 
MJ fragment) dominate as T increases. These trends 
reflect the more pronounced (quenching) effect of the 
electronic entropy on the shell-stabilization of the larger 
magic sizes (e.g., McJ" and M^j^). 

In the experiment [13], a distribution of relatively large 
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FIG. 1. Fission channels recorded in TOF measurements 
for Li^+, Na^+, and clusters, with iV = 24, 26, 28, and 
30. In each case, the peak on the right (marked by an arrow) 
corresponds to the parent cluster and the ones on the left of it 
correspond to the larger fission fragments. Fragment labelings 
are determined from the voltage setup. Peak heights are in 
arbitrary units. The splitting of the TOF ion-fragment peaks 
are due to the kinetic energy released during fission, which is 
of the order of 1 eV (see, e.g., the case = 26). 

neutral clusters is formed by a gas aggregation source. 
The clusters are ionized, photo-excited and warmed by 
a 15 ns laser pulse (Nd-YAG laser at hv = 3.50 eV). A 
rapid sequential evaporation follows the warming of the 
clusters for about 1 /xs (i.e., the residence time in the ion- 
izing region). This results in a charged cluster distribu- 
tion shifted down to smaller sizes, called an "evaporative 
ensemble" (EE) [13]. The cluster temperature associated 
with the EE is determined [13] by the cluster dissociation 
energy and by the experimental time window, i.e., the 
residence time in the ionizing/heating region. Under our 
experimental conditions, the temperatures of the EE's of 
lithium, sodium and potassium clusters are 700±100 K, 
400±100 K and 300±100 K, respectively. 

The EE of ionized clusters enters a tandem time-of- 
flight (TOF) mass spectrometer. Parents of interest are 
mass selected in the first TOF region. Because they 
are at a given temperature, subsequent unimolecular dis- 
sociation of these parents takes place and the resulting 
fragments are mass analyzed by the second TOF device. 
Enhanced optimization of the electrostatic potentials in 



the current experiments has led to significant increases 
in both the mass and energy resolutions. In particular, 
the latter allows us to resolve the kinetic energy released 
during the fission process (see the splitting of the peak 
of some of the fission fragments in Fig. 1). 

For cluster sizes in the vicinity of the appearence 
size (see below), the unimolecular dissociation of dou- 
bly charged clusters portrays the competition between 
the fission (la) and evaporation (lb) processes: 

M^+-^M+-|-M+_p (la) 
^ M|+ + Mn-p' , P' = 1 or 2 . (lb) 

We focus here on the fission decay channels. The ex- 
plored size domain is bounded by two values, iV_ and 
N^. The lower value iV_ is determined by the appear- 
ance size AT^PP (A^^PP = 24 for Li^+ and Na^ clusters, 
and 19 for clusters [4]), below which fission is the 
dominant dissociation process [see Eq. (la)]. Below this 
size no doubly charged clusters are present in the mass 
spectra, since their fast fission process prevents their ob- 
servation during the experimental time window. A_|_ is 
determined by the loss of the signal corresponding to the 
fission process, resulting from the increase of the inner 
part of the fission barrier with the size of the cluster rel- 
ative to the essentially constant monomer evaporation 
energy. Fig. 1 summarizes the results obtained for the 
three alkali metals in the same size range. 

Although we have recorded the fragmentation spectra 
for all the sizes 24 < A^ < 30, only results for clusters 
with an even number of atoms are shown (due to pairing 
they exhibit simpler features). In our experimental setup, 
only the heavy fragments Mp with P > N/2 are detected. 
The main features are as follows: For the lowest parent 
size (A^ — 24) and for all three elements, only one fission 
channel is present, i.e., the doubly magic + A/21. 
Increasing the parent size, however, results in differenti- 
ation between the three elements regarding the observed 
fission channels. A larger parent size (i.e., A^ > 26) opens 
the channels of heavier fragments and this trend is more 
pronounced for potassium than for lithium. In particu- 
lar, the trimer appears to be the dominant product in the 
Li fission, while Mj^ is the dominant one for K, with Na 
exhibiting an intermediate behavior. That is, the asym- 
metry of the fission process is reduced for the K clusters 
as compared to the Na and Li ones. 

We have performed FT-SCM calculations of the Q^- 
values for the various fission channels for all three alkali- 
metal sequences at the corresponding experimental tem- 
peratures. Since for each species the Coulomb repulsion 
(QR) at the scission point (measured by the energy re- 
lease) is approximately constant for all fission channels 
in the size-range considered here {X << 1), the effective 
inner part of the fission barrier can be approximated by 
adding the calculated Q-^-values to the QR term. Thus 
the ordering of the fission channels can be directly in- 
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FIG. 2. Theoretical SCM zero-temperature energetics 
(Q-values) for the processes Li^ — > Lip + Li^_p, A'^ — 24, 
26, 28, and 30. The dashed Unes correspond to LDM calcula- 
tions. In each case, the horizontal axis corresponds to various 
SCM channels 2 < P < iV - 2. 

ferred from the Q-^-values. 

The theoretical zero-temperature Q-values provide an 
inadequate description of the experimental data. For ex- 
ample, the Q-values for Li^ (plotted in Fig. 2) cannot 
explain the observed trend (see Fig. 1, left column) of 
strong asymmetric fission, with the Li^ being the pre- 
ferred channel for all parents iV = 24, 26, 28, and 30 [14]. 
In particular, the T = theoretically deduced preferred 
fission channels for = 26 and = 30 are (5,21) and 
(9,21), respectively, in contradiction with the measure- 
ments (see Fig. 1). In contrast, the theoretical Q"^-values 
(Fig. 3) agree well with the experiment. Specifically, for 
Lithium the calculations at T = 700 K show that Li^ 
is indeed the preferred channel for all parents; only for 
Ligfj", the Lijj^ channel becomes competitive, remaining 
however less abundant than the one involving the trimer. 

We reiterate that in Fig. 3 there is a clear tendency 
toward opening the channels involving the heavier magic 
fragments {P = 9 and P = 21), both when going up 
vertically from the lighter to the heavier parentss (from 

= 24 to = 30), as well as when going horizontally 
from left to right (i.e., from Li to K, with Na representing 
an intermediate case). 

The physics underlying this agreement between the ex- 
periment and the FT-SCM results can be revealed by an 
examination of the influence of the electronic entropy on 
the various fission channels {P, N ~ P) as a function of 
P. For brevity, we limit our discussion to the case of the 
Kg^ parent; the conclusions, however, extend to all three 
species and the other parent sizes. Figure 4 displays the 
FT-SCM Q^-values for at five temperatures (200 K 
to 1000 K). At the lowest temperature, the (9,21) chan- 
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FIG. 3. Finite-temperature SCM results (Q^^-values) for 
the fission process M^ Mp + M'^_p for Li^, Na^, and 
with N = 24, 26, 28, and 30. The corresponding tem- 
peratures (in degrees Kelvin), taken from the experiments, 
for each species are shown at the top. The dashed lines corre- 
spond to LDM calculations. In each case, the horizontal axis 
corresponds to various SCM channels 2 < P < N — 2. The 
length of the vertical bar in each frame corresponds to 0.4 eV. 

nel is the preferred one; however, an increase in T is ac- 
companied by quenching of the shell effects; i.e., at high 
temperature (1000 K) the shell effects have almost van- 
ished, and the Q^ curve exhibits a smooth shape similar 
to that of the LDM (although with a steeper variation) 

Most important for the competition between the var- 
ious magic fragments, however, is the fact that the 
electronic-entropy quenching of the shell effects is non- 
uniform, with the higher occupied shells of the heavier 
fragments being influenced the most as T increases, since 
the energy spacing between the major electronic shells, 
hiUoiN) = 49 eV • a'^/{r'^^N^/^), decreases with larger A^. 
Consequently, a higher T favors the trimer {M'^) channel 
over those involving the magic Mg and fragments. 
This happens even in the case of A^ = 30 where the com- 
peting (9, 21) channel is doubly magic. Note that, unlike 
the T = 700 K experimental temperature of Li, the tem- 
peratures of 300 K for Potassium and 400 K for Sodium 
are low enough so that the heavier-fragment channels re- 
tain their T = original prevalence for all three parent 
sizes A^ = 26, 28, and 30. 
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FIG. 4. Theoretical FT-SCM temperature dependence of 
the Q"^-values for the process Kg J Kp + K3'Q_p. The 
dashed hnes correspond to LDM calculations. In each case, 
the horizontal axis corresponds to various SCM channels 
2 < P < 28. In the bottom right frame, we display fsh vs. r 
(see text) for the Kg" and Kg" fragments. 



Further insight can be gained by examining the ra- 
tio U{N) = AF3h(7V,T)/|AF,h(Af,0)| of the free-energy 
shell-correction term at r = 2Ti'^kBT /hu)Q{N) over the 
shell correction term at r = 0. For a schematic single- 
particle spectrum consisting of equally spaced major 
shells with the same degeneracy, it can be shown in the 
case of closed shells (see p. 608 of Ref. [2] and Ref. [10]) 
that fsh{N) ~ — rexp(— r) for r > 1. Thus, because 
of the large prefactor 27r^, the electronic-entropy effects 
manifest themselves even at the rather low temperatures 
of the cluster experiments. In Fig. 4 (bottom frame on 
the right), we compare the /sh's for and calcu- 
lated with the realistic single-particle spectrum of our 
FT-SCM. Note that, as a function of T, the realistic spec- 
trum suppresses the shell effect of the heavier-than-the- 
trimer magic fragments at a faster rate than the estimate 
above. 

The controlling influence of the electronic entropy on 
the stabilization of electronic shells, discussed previously 
[6,10] in studies of ground-state properties of clusters, 
has been shown here to govern the fission patterns and 
their temperature dependence at the low-fissility regime 
[15]. This provides the impetus for further temperature- 
controlled studies of charge instabilities in clusters. 
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